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Abstract. Event shapes and jet shapes in neutral current deep inelastic scattering and photoproduction are 
studied with the HI and ZEUS detectors at HERA. The measurements are compared to next-to-leading- 
order QCD calculations and Monte Carlo models. The strong coupling constant is determined from subjet 
multiplicities. 



1 Introduction 

The HERA collider is a precision tool for the study of 
QCD. The HI and ZEUS Collaborations have performed 
a wealth of analyses which investigate the hadronic final 
state in ep collisions. In addition to the measurement of 
plain jet cross sections, it is rewarding to study the sub- 
structure of jets, or in the case of event shapes, the struc- 
ture of the whole hadronic final state. In this case, uncer- 
tainties connected to the determination of cross sections 
for jet production do not compromise the results. 
In the regime of high transverse momenta, experimental 
uncertainties and non-perturbative effects are small, which 
makes possible precision tests of perturbative QCD and 
the extraction of the strong coupling constant. 
On the other hand, understanding the nature of hadroni- 
sation is a challenge, which is important for all high energy 
experiments involving hadrons. In this area, the concept 
of power corrections offers prospects beyond phenomeno- 
logical models implemented in Monte Carlo generators. 



where fig b j(y C ut) is t ne number of subjets in jet i and iVj e t s 
is the total number of jets in the sample. 
The ZEUS collaboration has measured the mean subjet 
multiplicity in neutral current deep inelastic scattering 
(NC DIS) as a function of E^ E3- Jets were defined by 
using the longitudinally invariant ki cluster algorithm in 
the laboratory frame. In Fig.rjthe measurements of (ro s bi) 
are compared with the NLO prediction from DISENT 2 , 
corrected for hadronisation cfccts. The overall description 
of the data is good, so that the measurements can be used 
to make a determination of a s . The region used for the 
fit was restricted to E 1 ^ > 25GeV, where hadronisation 
corrections differ by less than 17% from unity. 
The QCD fit of the mean subjet multiplicity for 25 < 
El? < 71 GeV at y cut = 1CT 2 yields 
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This result is consistent with recent determinations by the 
HI jHEl and ZEUS |5K| Collaborations and with the PDG 
value [7j- The total error is dominated by the theoretical 
uncertainty due to missing terms beyond NLO. 



2 Subjet Multiplicities 

The internal structure of jets is determined by QCD ra- 
diation processes. Non-perturbative fragmentation effects 
become small at high transverse jet energies Et and values 
of the resolution scale y cn t that are not too low. Pertur- 
bative QCD calculations can be applied. In this region of 
phasespace they depend only weakly on parton density 
functions, but are directly proportional to the strong cou- 
pling constant a s . 

The mean subjet multiplicity is defined as 

(n S bj(2/cut)) = TT Y] «sbj(2/cut), 



3 Substructure Dependence of Jet Cross 
Sections 

Quarks and gluons carry different colour charges, there- 
fore the internal structure of a jet is expected to depend 
on the type of the primary parton from which it origi- 
nated. Obervables which probe the internal structure are 
the subjet multiplicity (section |2J) and the integrated jet 
shape &(r), which is defined as 



E T (r) 
E^ ' 



On behalf of the HI and ZEUS Collaborations 



where Et(t) is the transverse energy within a cone of ra- 
dius r — \J A(\> 2 + Ar/ 2 about the jet axis. 
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Fig. 1. The mean subjet multiplicity of an inclusive jet sample 
in NC DIS, corrected to hadron level, as a function of Erjet- 



In photon-proton reactions, two types of processes con- 
tribute to jet production. In the direct process the photon 
interacts with a parton in the proton as a whole, whereas 
in the resolved process it behaves as a source of partons 
which subsequently interact with partons from the pro- 
ton. Jets from direct and resolved processes exhibit differ- 
ent angular distributions and Et spectra. Consequently, 
dijets from both processes differ in their distributions of 
invariant mass M 33 and the angle between the jet-jet axis 
and the beam direction in the dijet centre-of-mass system 
cos (9*. The cross section for subprocesses with gluons in 
the final state is larger for the resolved than for the direct 
process, if an appropriate region of phase space is chosen. 
Therefore, the experimental investigation of jet substruc- 
ture as a function of jet kinematics allows one to study 
the characteristics of quark- and gluon-initiated jets. 
The ZEUS Collaboration has measured differential inclu- 
sive jet and dijet cross sections in photoproduction [SJ. 
The internal structure of jets, in terms of the integrated 
jet shape and subjet multiplicity, has been used to select 
samples enriched in quark- and gluon-initiated jets, clas- 
sified as "thin" and "thick" samples, respectively. Fig. [21 
shows the measured inclusive jet cross section as a func- 
tion of Et and 77j 0t , and the dijet cross section as a func- 
tion of | cos 9* \ and M 33 . The distributions are in general 
well described by calculations using PYTHIA, which have 
been normalised to the total measured cross section. The 
"thick" data sample behaves according to expectations for 
gluon initiated jets in the ivestigated phase space: these 
jets are more forward, have a softer E t spectrum (inclu- 
sive jets), lower invariant mass and the jet axis is closer to 
the beam direction (dijets). This behaviour is consistent 
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Fig. 2. Differential jet cross section in photoproduction for 
an inclusive (upper) and dijet (lower) sample. The jets were 
assigned to "thick" and "thin" subsamples according to their 
shape and subjet multiplicity. 



with the expected dominance of the resolved subprocess 
for gluon initiated jets. 



4 Event Shapes 

Event shape variables are designed to study QCD by mea- 
suring properties of the flow of hadronic energy-momen- 
tum. A suitable frame of reference in DIS is the Breit 
frame, which divides an event into hemispheres correspond- 
ing to the proton remnant and the hadronic final state 
evolving from the struck parton. Event shape observables 
are then usually defined for the particles in the current 
hemisphere (CH) alone, thereby avoiding non-perturbative 
effects connected with the proton remnant. A commonly 
used event shape variable is thrust r = 1 — T: 



T 



E/ l£ CH \Ph\ 



where ph is the momentum of particle h in the current 
hemisphere and p z h this momentum projected on the axis 
of the exchanged virtual photon. The description of such 
infrared-collinear safe observables by fixed order calcula- 
tions faces two difficulties. Firstly, non-perturbative hadro- 
nisation effects can be large, even at formally perturbative 
scales. In the approach by Dokshitzer, Webber et al. [H] 
these effects can be addressed using power corrections pro- 
portional to l/Q, governed by one universal free parame- 
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ter an. Secondly, the convergence of the perturbative se- 
ries at low values of the event shape variables is very poor, 
making a resummation to all orders in a s necessary. Re- 
summed NLL calculations, matched to NLO, have recently 
become available for two kinds of thrust, r c and r, the jet 
broadening B, the jet mass po and the C-parameter [TH] . 
While these event shapes are mainly sensitive to the tran- 
sition from 1-jet like to 2-jet like topologies, there are now 
also measurements and calculations available for the out- 
of-plane momentum K out and the azimuthal correlation 
X, which probe the transition from 2-jet to 3-jet topolo- 
gies [TT] . 

Mean values of event shapes as a function of Q have been 
studied by the HI [H] and ZEUS [T3] Collaborations. The 
HI Collaboration has also measured spectra of event shape 
variables :14 . Fig.[31shows the distributions of thrust and 
out-of-plane momentum in bins of momentum transfer Q. 
The thrust spectra are compared with the results of a 
fit based on resummed calculations matched to NLO and 
including power corrections. For the out-of-plane momen- 
tum the matching is not yet available, therefore a com- 
parison with the RAPGAP Monte Carlo event generator 
1 15) is presented. This includes LO matrix elements and 
parton showers. Both observables are well described over 
the whole region of phase space. 

The ZEUS Collaboration has made a fit to the mean values 
of event shape variables using DISASTER++ ^H] and the 
Dokshitzer- Webber power corrections, as shown in Fig. [J] 
(left). The contours correspond to 1-cr statistical errors. 
Consistent values of a s are obtained for four of the shape 
variables with an values that agree within 10%. However, 
the lack of overall consistency in a s and an suggests that 
higher order processes which are not included may be im- 
portant. On the right hand side of Fig. HI results °f a nt to 
event shape spectra by the HI Collaboration are shown, 
including resummed calculations. Here the contours in- 
dicate 1-cr statistical and experimental systematic errors, 
added in quadrature. The results for the five studied event 
shapes are consistent with each other and with the world 
average of a s |17| . The universal non-perturbative param- 
eter an is confirmed to be 0.5 at the 10% level. 



5 Conclusion 

Differential distributions of jet and event shapes have been 
measured by the HI and ZEUS Collaborations with high 
precision. Calculations based on perturbative and non- 
perturbative QCD describe the measurements over a large 
region of phasespace. 
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Fig. 3. Measured distributions of 1-thrust r (left) and out-of- 
plane momentum K out (right), r is compared with the results 
of a fit based on NLO QCD including resummation and power 
corrections, K ou t is compared with results from the RAPGAP 
Monte Carlo generator. 
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Fig. 4. 1-ct contours in the (a 3 ,ao) plane from fits to event 
shape mean values (left) and differential distributions (right). 
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